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Abstract. We studied by X-ray reflectivity the behaviour of fully hydrated solid-supported lipid multilayers 
under the influence of a transverse electric field, under conditions routinely used in the electroformation 
process. The kinetics of sample loss (unbinding) was measured as a function of the amplitude and frequency 
of the applied field by monitoring the integrated intensity of the Bragg peaks. We also performed a time- 
resolved analysis of the intensity of the first Bragg peak and characterized the final state of the sample. 

PACS. 61.10.Kw X-ray reflectometry (surfaces, interfaces, films) - 87.16.Dg Membranes, bilayers, and 
vesicles - 87.50.Rr Effects of electric fields on biomolecules, cells and higher organisms 


1 Introduction 

The influence of an external electric field upon lipid bilay¬ 
ers is of paramount importance, both from a practical and 
a fundamental point of view. Cell electroporation [Tj has 
been long employed for introducing exogenous molecules 
into cells, while electroformation 3| is the best method 
to date for obtaining giant unilamellar vesicles, which can 
be used as model systems for the cell membrane and as 
transfection factors (see IJj and references therein). Al¬ 
though very important, these phenomena are still in need 
of a clear-cut explanation. This lack can be understood in 
light of the fact that both processes are instabilities involv¬ 
ing large values of the electric field. As such, they lead to 
large deformations, irreversibly changing the system struc¬ 
ture : creating pores in the membrane (electroporation) or 
‘peeling’ the bilayers off the substrate (electroformation). 

Here, we investigate the electroformation process [3], 
which involves applying an AC electric field across a lipid 
film deposited on a substrate. The later stages of the 
process, namely the formation of small vesicles and the 
’’ripening” leading to the appearance of giant vesicles have 
already been studied in detail [3|. However, the very first 
step, during which the lipid bilayers unbind from the sub¬ 
strate, has not been quantified. In the present work, we 
focus on characterizing this initial step, using the X-ray 
reflectivity technique [5] , which has the advantage of only 
being sensitive to the bound membranes [6]. Thus, we can 
determine the quantity of lipid still left on the substrate at 
a certain time (averaged over the fotprint of the beam) as 
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well as structural changes at the microscopic level (lamel¬ 
lar spacing, fluctuation amplitude etc.) that are inaccessi¬ 
ble by light microscopy. It should be noted that the recent 
experimental work of Burgess et al. [7: studies the influ¬ 
ence of the (static) surface charge on the structure of a 
unique supported bilayer, using neutron reflectivity. 

Several theories were put forward to explain different 
aspects of the electric field effect upon the stability of lipid 
membranes (flat [5! or under the form of vesicles [9lflT)|. 
charged mmm or uncharged ED- To our knowledge, 
the only theoretical paper relevant to our experimental 
configuration (uncharged flat membrane under the influ¬ 
ence of a transverse electric field) is the one by Sens and 
Isambert m- They predict an undulation instability of 
the membrane, with a typical wavelength in the submi¬ 
cron range. 


2 Experimental method 

2.1 Sample preparation and environment 

The field unbinding experiments were performed using the 
zwitterionic lipid l,2-dimyristoyl-sn-glycero-3-phosphocho- 
line (abbreviated in the following to DMPC), bought from 
Avanti (Alabaster, AL, USA) and used without further 
purification. The lipid was dissolved in a 1:1 (vol/vol) 
mixture of chloroform and TFE (2- 2-2-trifluoroethanol) 
at concentrations between 5 and 20 mg/ml. An amount 
of 0.1-0.2 ml of the solution was pipetted onto carefully 
cleaned silicon substrates of a size of 15 x 25 mm 2 cut 
from conductive (highly doped) commercial silicon wafers 
(Silchem Gmbh, Freiberg, Germany). Depending on the 
desired film thickness, the samples were either left to dry 
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slowly, in which case all the lipid remains on the sub¬ 
strate, yielding stacks of 1000-3000 bilayers (depending 
on the solution concentration and volume), or submitted 
to the spin-coating process m- In this case, immediately 
after deposition the substrate was accelerated to rotation 
(3000 rpm), using a spin-coater. In this case, the number 
of bilayers 5 < N < 25 can be controlled. In both cases, 
the samples were left to dry under a laminar flux hood 
at room temperature for a few hours and then exposed to 
high vacuum at 40 °C overnight to remove any remaining 
traces of solvent. They were finally stored at 4 °C until the 
measurement. 

The experiments were performed in a plexiglas cham¬ 
ber with kapton windows (Fig. [l]). The sample contacts 
and the counter electrode are in stainless steel. The elec¬ 
tric field is applied using a function generator (Agilent 
Technologies). The peak-to-peak amplitude is given in the 
text as Vpp. The cell was mounted on a metal heating 
stage, temperature-controlled by water flow from a heat¬ 
ing bath (Julabo Gmbh, Seelbach, Germany). After mount¬ 
ing the sample, the hydrating solution was gently added 
to avoid washing the lipid film off the substrate. We used 
ultrapure water (Millipore, Bedford, MA). 



The data was corrected for decreasing electron ring cur¬ 
rent and for the diffuse contribution (by subtraction of an 
offset scan). Finally, an illumination correction was per¬ 
formed. 


2.3 Time correlation function 

For some samples we determined the autocorrelation func¬ 
tion of the scattering signal. After aligning the sample 
(usually with the first Bragg peak in specular condition), 
the signal from the scintillation detector was fed into a 
multiple-tau correlator (ALV-5000/EPP, from ALV Gmbh, 
Langen, Germany) and the corelation function [IT] g(t ) = 

— ' l/r w as obtained. g(t ) was averaged over 300 s. 
We performed these measurements both on thick (spread) 
and thin (spin-coated) samples. No detectable modulation 
was present for thick samples; however, for thin samples 
and at moderate frequencies, g(t ) exhibited an oscillation 
at the frequency of the applied field. When this feature was 
observed, care was taken to rule out any artefacts : we 
checked that the modulation disappeared when unplug¬ 
ging the cable and when the sample chamber was vigor¬ 
ously flushed with solvent to remove the sample from the 
substrate, maintaining the experimental configuration and 
the applied field. This latter test conclusively shows that 
the oscillation is indeed a feature of the lipid bilayers, and 
not of the substrate or mounting. 


3 Results 
3.1 Field amplitude 

We studied the unbinding kinetics as a function of the 
amplitude of the applied electric field. All the measure¬ 
ments discussed in this subsection were performed at a 
fixed field frequency v = 10 Hz and at a temperature 
of 35 °C. First, we present results obtained at low am¬ 
plitudes, U < 2 Vpp, on spin-coated samples. Figure [2] 
shows a typical reflectivity spectrum before applying the 
held (A). After the held is turned on, only the vicinity 
of a Bragg peak is scanned, as shown in the inset (B). 
For all the curves, the measurements were performed on 
The X-ray reflectivity measurements presented here were the second Bragg peak, where the background (due to 
carried out at the bending magnet beamline D4 of HASYLABfehe interference with the substrate) is hatter and can be 

more easily subtracted. The area under the Bragg peak is 
then integrated and corrected for the time variation of the 
primary beam. Since the integrated intensity is roughly 
proportional to the sample volume (number of layers x 
illuminated area), we can thus monitor the kinetics of un¬ 
binding averaged over a few mm 2 of sample area. The 
integrated intensity is shown in panel (C) as a function of 
time, for different values of the applied field U = 0.2, 0.3, 
0.45, 1.0 and 2.0 Vpp. The longest evolution measured (at 
0.2 Vpp) was recorded over almost ten hours. 

For the lowest applied fields, the intensity decay with 
time has a sigmoidal shape, characterized by a ’’waiting 
time” and a ” decay time”; we quantified this tendency by 
fitting the curves with a hyperbolic tangent : 


2.2 X-ray reflectivity 


DESY in Hamburg, Germany. At D4, a single-reflection 
Si(lll) monochromator was used to select a photon en¬ 
ergy of 20 keV, after passing a Rh mirror to suppress 
higher harmonics. Primary beam intensity is around 5 10 8 
counts/s. The use of high energy X-rays has two reasons : 
it reduces the beam damage to the sample and minimizes 
absorbtion in the solvent (the transmission factor is about 
0.25 for 2 cm of water). 

The chamber was mounted in horizontal scattering ge¬ 
ometry, and the reflected beam was measured by a fast 
scintillation counter (Cyberstar, Oxford), using computer- 
controlled aluminum absorbers which attenuate the beam 
at small q z to prevent detector saturation. Incident and 
exit beams were defined by a system of motorized slits. 
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Fig. 2. A : X-ray reflectivity curve of a thin DMPC sample on Si substrate. B : Zoom of the region around the second Bragg 
peak (dashed line in A), at different times, under an applied electric field of 0.45 Vpp and 10 Hz. C : Integrated intensity of the 
second Bragg peak as a function of time, for different amplitudes of the electric field. For the three lowest voltages, the fit with 
equation © is also shown (solid line). 


m = 


1 — tanh 


t-t 0 

At 


yielding the following values : 


(1) 


u [Vpp] 

to [s] 

At [s] 

Sample 


0.2 

20000 ± 300 

8000 ± 500 

thin 

Fig.EJC 

0.3 

10000 ± 40 

5500 ± 60 


55 

0.45 

2500 ± 70 

1800 ± 120 

» 

55 

1.0 

0 [fixed] 

1500 ± 250 

II 

55 

2.0 

0 [fixed] 

800 ± 100 

» 

55 

5.0 

0 [fixed] 

1300 ± 100 


not shown 

1.0 

3200 ± 80 

1500 ± 120 

thick 

not shown 


Table 1. Parameters of the intensity decay as a function of the 
amplitude of the field (at a fixed frequency u = 10 Hz) fitted 
with equation© 


At higher field amplitudes (above lVpp), the ’’wait¬ 
ing time” vanishes, so that we set tp = 0. The remaining 
decay could be equally well described by an exponential 
function. 

We also studied thick samples under the same experi¬ 
mental conditions. For comparison, the evolution of the 
scattered intensity was fitted using the same model as 
for the thin samples (equation [0 and the parameters are 
given in Table 1. 

The main interest of using thick samples is, however, 
the possibility of studying the diffuse scattering signal and 
of obtaining information about the bilayer position fluc¬ 
tuations. Extensive measurements were performed, both 
at D4 (HASYLAB) and at the ID1 beamline at the ESRF 
(Grenoble, France); no difference was detected in the dif¬ 


fuse spectrum between the situation with or without elec¬ 
tric field (data not shown). We therefore conclude that the 
bilayers within the sample are not significantly affected by 
the presence of the field (no enhancement in the fluctua¬ 
tion amplitude). 


3.2 Comparison with thermal effects 

The dependence of the unbinding rate on the amplitude 
of the electric field is strikingly similar to the kinetics 
recorded as a function of the temperature 118l . The same 
pattern of a plateau followed by a decay (at lower tempera¬ 
ture) and the disappearance of the plateau (at higher tem¬ 
peratures) emerge, albeit on a completely different time 
scale (see Figure©. 

The lipid used was l-oleoyl-2-palmitoyl-sn-glycero-3- 
phosphocholinc (OPPC), bought from Avanti Lipids, Al¬ 
abama. The samples were prepared as for the electric field 
studies (by spreading), except for using isopropanol as a 
solvent and for applying a hydrophilizing treatment to the 
substrates, by washing them in a saturated KOH solution 
in ethanol for about a minute. Subsequently, they were 
rinsed several times with ultra-pure water. 

The data has been taken in a temperature-controlled 
steel chamber sample immersed in excess ultra-pure water 
(the beam path in water is 15 mm), on an in-house setup 
using a sealed tube generator (Mo K a radiation, E = 
17.4 keV) and a two-circle goniometer (Siemens D500). 
The integrated diffuse scattering of the first Bragg sheet 
was used as a measure of the total scattering volume. 

No high-temperature measurements were performed 
for DMPC, so that a direct comparison with the results of 
subsection I3J1 is not possible. However, several measure¬ 
ments were made using one DMPC sample, at tempera- 
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Fig. 3. Integrated intensity of the first Bragg sheet (diffuse 
signal) of an OPPC sample as a function of time, for different 
temperatures. For the three lowest temperatures, the fit with 
equation © is also shown (solid line). 


tures between 22 and 30 °C, over a period of four weeks, 
without any apparent sample loss. 




3.3 Field frequency 

The temperature was set at 30 °C for these measurements. 

A representative example is shown in Figure © graph 
A shows the specular scan (corrected for the offset contri¬ 
bution) and the adjustment with the model given in [191 
ED] (red solid curve), as well as the associated Fresnel re¬ 
flectivity (dashed line). The inset shows the same data di¬ 
vided by the Fresnel reflectivity. From the model, we infer 
that the sample consists of four layers, with approximate 
coverage ratios (starting at the substrate) : 1.0, 0.8, 0.45, 
0.15 (other sample parameters, such as the amplitudes of 
the different Fourier components of the electron density 

7T k]3 -7" 1 

profile and the Caille parameter g = — ———— quantifying 

_ 2 BXd z 

the amplitude of the thermal fluctuations were taken from 
previous measurements under the same conditions [20]). 
Graph B shows a time correlation function measured on 
top of the first Bragg peak ( q z = O.lOlA , correspond¬ 
ing to a repeat spacing d = 62. 0A), under an applied 
electric field with an amplitude of 5 Vpp and a frequency 
v = 30 Hz. The correlation function is well described by a 
decaying oscillation at the applied frequency : 

g{t) — 1 = go cos(27ri4) exp(— t/r) (2) 

where the amplitude go ~ 5 10~ 3 and the decay time r = 
0.23 s. Similar measurements were performed at slightly 
different frequencies; the behaviour of g(t) follows the ap¬ 
plied frequency, ruling out the presence of an external os¬ 
cillator. As discussed above (subsection 12.3|h we took ad¬ 
ditional step to confirm that the signal is indeed a genuine 
contribution of the sample. 

We recorded the correlation function on top of the first 
Bragg peak for increasing frequencies and noticed that, 


Fig. 4. A : Reflectivity of the spin-coated sample discussed in 
the text (open symbols) adjusted with the model described in 
[191120] (red solid curve) and the associated Fresnel reflectivity 
(dashed line). The inset shows the same data scaled by the 
Fresnel reflectivity. B : Correlation function measured on the 
first Bragg peak ( q z = O.lOlA ) under an applied electric 
field with v = 30Hz (open symbols) and fit with a damped 
cosine (red solid curve). See text for details. 


above a threshold of the order of i/ max = 100 Hz, estimated 
from measurements performed on several samples, g(t) no 
longer exhibits oscillations at the applied frequency. 

In the range where the correlation function can be 
properly measured (about 10 to 50 Hz), the recorded fre¬ 
quency v is very precisely that of the applied electric field. 
The decay time exhibits an intriguing dependence on v : 
r ~ 7/zq but more systematic measurements would be 
needed to confirm whether this is an intrinsic feature of 
the sample or an artifact of the setup. 

The frequency of the applied field also has a strong 
effect on the unbinding kinetics, as can be seen in Figure 
[5] The measurement was measured on a thick (about 3000 
bilayers) sample. We started by applying a 10 Vpp field 
at 1 kHz and then decreased the frequency at 100 Hz and 
finally at 10 Hz. The change in the decay rate of the peak 
is very strong. In addition, an abrupt decrease in intensity 
was observed when switching from 100 to 10 Hz; this effect 
is systematically observed in all samples studied at high 
field amplitude (data not shown), and could be related to 
the unbinding of weakly adhering surface patches. 
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Fig. 5. Unbinding kinetics for a thick (spread) sample, as a 
function of the frequency of the applied electric field. Note the 
logarithmic intensity scale. 



Fig. 6. Reflectivity of a spin-coated sample in the final state, 
after applying a strong electric field (open symbols) adjusted 
with the model described in [20] (red solid curve) and the as¬ 
sociated Fresnel reflectivity (dashed line). The inset shows the 
same data scaled by the Fresnel reflectivity. 


3.4 Final state 

An interesting point to address is also the final state, de¬ 
fined as the point at which applying a reasonably strong 
electric field no longer has any observable effect on the 
sample. In most of the measurements performed on thin 
(spin-coated) samples, it appears that only one or two bi¬ 
layers remain on the substrate. An example is shown in 
Figure [G] where a sample prepared according to the pro¬ 
tocol in m using a DMPC concentration of 2 mg/ml was 
exposed to an AC field with an amplitude of 10 Vpp at a 
frequency of 10 Hz, and afterwards at 5 Hz for more than 
10 minutes in each case; the final reflectivity (open sym¬ 
bols) can be described by assuming that only two bilayers 
are left on the substrate, with coverage ratios of 0.8 and 


0.3, starting at the substrate. This spectrum being a rep¬ 
resentative one for spin-coated samples, we can conclude 
that the bottom layer (closest to the substrate) is largely 
unaffected by the electric field. 

4 Discussion and Conclusion 

At the beginning of this study, two kinds of results were 
expected : long-time effects (unbinding kinetics under elec¬ 
tric field) and changes in the short-time behaviour (the 
fluctuation spectrum) of the lamellar phase. 

Unbinding kinetics : The unbinding of lipid bilay¬ 
ers from a solid substrate represents the first step of the 
electroformation process. We determined the kinetics of 
unbinding under field for various applied voltages, both 
for thin (tens of bilayers) and thick (thousands of bilay¬ 
ers) samples. A two-step process was observed at low volt¬ 
age, consisting in a very slow decay (’’incubation” plateau) 
followed by a fast sample loss. At higher applied voltages, 
the first step disappears. Surprising similarities with the 
previously characterized thermal unbinding kinetics were 
thus revealed. The frequency of the applied field is also 
important: the unbinding no longer occurs when the fre¬ 
quency exceeds a few hundred Hz. 

The field-induced instability discussed here and the 
disintegration pathway of the lamellar stack has been cor¬ 
roborated in many sample series of DMPC, at varied field 
amplitudes, and frequencies, both for thin and thick stacks 
m- Moreover, it was found that the presence of an os¬ 
motic stressor (polyethylene glycol, MW 20000) in solu¬ 
tion does not prevent the instability, but does slow down 
the decay rate 121) ■ However, the results of measurements 
on thick stacks and/or in the presence of a osmotic stress 
being rather noisy, we do not discuss them here. 

Fluctuations : As we know that the applied field in¬ 
duces the unbinding, it must necessarily have an effect 
at the microscopic scale (i.e. that of the individual mem¬ 
brane). This effect could be at the level of the structural 
parameters (changes in the bilayer and/or water layer 
thickness) but could also involve enhanced fluctuations 
of the bilayers, as invoked in the case of thermal unbind¬ 
ing. In addition, recent theoretical work m has suggested 
that a transverse electric field could induce an undulation 
instability of a (low conductivity) membrane. 

We could detect no change in the structural parame¬ 
ters or in the diffuse scattering for thick samples; in thin 
samples, on the other hand, the intensity of the Bragg re¬ 
flection exhibits a modulation at the frequency of the ap¬ 
plied field. This leads us to the tentative conclusion that 
the electric field only affects a few layers at the top of 
the stack and that, in thick samples, their contribution to 
the scattered signal is swamped by that of the bulk. It 
appears that we are dealing with a surface phenomenon, 
which is reasonable in light of the fact that the electrofor¬ 
mation phenomenon yields unilamellar vesicles; at some 
point, the bilayers must peel off one by one. 

This conclusion is not in contradiction with the theo¬ 
retical model of Sens and Isambert m, since they also ob¬ 
serve a strong change in dynamics in the presence of fixed 
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boundaries. It would be reasonable to assume that motion 
of the membranes within the stack is hindered by their 
neighbours (and the boundary condition at the sustrate), 
while the top ones are relatively free to fluctuate. A more 
quantitative comparison cannot be made, since we cannot 
determine accurately the number of affected bilayers at 
the surface, the fluctuation amplitude and length scales; 
it should also be noted that membrane integrity is a key 
ingredient of their model (the electric current flows across 
the membrane, and the transverse electric field is ampli¬ 
fied in proportion to the resistivity of the membrane). 

This condition is not fulfilled in our samples, where 
defects (dewetted patches) are known to appear [22]. The 
electric current is more likely to flow through the defects, 
rather than across the bilayers; this is also the case of elec¬ 
troformation procedures, where the deposited lipid film 
need not cover the entire surface of the electrode [3] . 

We conclude that the presence of the electrodynamical 
instability studied by Sens and Isambert [15l is neither 
confirmed nor can it be refuted by our experiments; a 
more careful setup, involving e. g. a freely suspended lipid 
film, would be needed. 

In these conditions, an alternative explanation could 
be related to the presence of electro-osmotic flow in the 
solution, but more comprehensive tests must be performed 
in order to test this hypothesis. 
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